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Abstract 

We investigate the consequences of enforcing local color neutrality on the color superconducting 
phases of quark matter by utilizing the Nambu-Jona-Lasinio model supplemented by diquark and 



0^ \ the t'Hooft six-fermion interactions. In neutrino free matter at zero temperature, color neutrality 

guarantees that the number densities of u, d, and s quarks in the Color-Flavor-Locked (CFL) 



phase will be equal even with physical current quark masses. Electric charge neutrality follows as a 
consequence and without the presence of electrons. In contrast, electric charge neutrality in the less 



o 

(N 

' symmetric 2-flavor superconducting (2SC) phase with ud pairing requires more electrons than the 
(N ' 

• normal quark phase. The free energy density cost of enforcing color and electric charge neutrality 



in the CFL phase is lower than that in the 2SC phase, which favors the formation of the CFL 



Qh' phase. With increasing temperature and neutrino content, an unlocking transition occurs from the 

' CFL phase to the 2SC phase with the order of the transition depending on the temperature, the 



quark and lepton number chemical potentials. The astrophysical implications of this rich structure 
in the phase diagram, including estimates of the effects from Goldstone bosons in the CFL phase, 
are discussed. 
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Studies of QCD at high baryon density have led to the expectation that quark matter 
is a color superconductor in which the pairing gaps of unlike quarks [ud, us, and ds) are 
as large as 100 MeV. For three massless flavors, a symmetric ground state called the Color- 
Flavor-Locked (CFL) phase, in which BCS-like pairing involves all nine quarks, is favored 
||T|, 0. At lower density and for physically relevant values of the strange current quark mass 
(100 < ms/MeV < 300), a less symmetric (2SC) phase in which only the light up and down 
quarks {mu,d < 10 MeV) pair is expected P, ^. For recent reviews, see Refs. 

With the exception of the work by lida and Baym 0], and more recently by Alford and 
Rajagopal [0, little attention has been paid to the issue of color neutrality in superconducting 
quark phases. These works are the primary motivation for this study. The issues addressed 
in this work are similar to those addressed by Alford and Rajagopal|0 who perform a model 
independent analysis that is valid when <C fi and A ~ rn^/fJ', where A is the pairing 
gap and fi is the quark number chemical potential. We employ an extended version of the 
the Nambu-Jona-Lasinio model (NJL hereafter), which shares many symmetries with QCD 
including the spontaneous breaking of chiral symmetry, and calculate the thermodynamic 
potentials, Q, and pairing gaps. A, self-consistently in the CFL and 2SC phases. Our analysis 
leads to results that complement some of the conclusions in Ref. |^. There are, however, 
several aspects in which we go further. First, we employ a self-consistent model which 
uniquely determines both the diquark and the quark-anti quark condensates. Second, since 
the realization of color and electric charge neutrality becomes non-trivial only for physically 
relevant values of m^, we retain terms to all orders in nis in our calculation of Q. As noted 
in Ref. 10], this is particularly important for understanding the phase structure of quark 
matter at densities (or equivalently, /i) of relevance to neutron stars, since m^/ /i is not small 
compared to unity. In addition, we establish the phase structure of superconducting color- 
neutral quark matter at finite temperature and lepton content which was not considered in 
@, 1^, but is relevant for studies of proto- neutron stars. 

Charges, chemical potentials, and color neutrality: Bulk, homogeneous matter 
must be neutral with respect to charges which interact through the exchange of massless 
gauge bosons. Otherwise, the free energy density cost would be infinite. In the CFL phase, 
diquark condensation breaks color symmetry and all eight gluons become massive via the 
Higgs mechanism. Similarly, in the 2SC phase, SU{3)c is broken down to SU{2)c and five 
of the eight gluons become massive. In both the CFL and 2SC phases, however, an U{1) 
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gauge symmetry remains unbroken The associated charge is called Q. The generator for 
this charge in the CFL phase is a linear combination of the usual electric charge Q and a 
combination of color generators T3 and Tg, and is given by 

where Q = diag{2/3, —1/3, —1/3) in flavor space, and T3 = diag{l, —1,0) and Tg = 
diag{l/\/3,l/\^,—2/\/3) in color space. 

The color superconducting phase is, by construction, neutral with respect to Q charge. 
Why then should we impose, in addition, local color neutrality? As noted earlier, gluons 
become massive in the superconducting phase and the free energy density cost of realizing a 
non-zero color density in bulk matter need not be infinite. Further, although a finite sample 
embedded in the normal state must be a color singlet, this alone does not require local color 
neutrality since color singletness is a global constraint. Hence a heterogeneous phase with 
colored domains of typical size similar to the color Debye screening length is a possibility. 
However, in a homogeneous and color conducting medium a color charge density in the bulk is 
unstable as it generates a chromo-electric field resulting in the flow of color charges 0. Color 
neutrality is therefore a requirement for the homogeneous phase. Neutrality with respect to 
charges associated with T3 and Tg is achieved by introducing appropriate chemical potentials 
/is and yUg in analogy with the charge chemical potential /xg. As noted in Refs. [0, ||, color 
neutrality is a prerequisite for color singletness, but the additional free energy density cost 
involved in projecting out the color singlet state is negligible for large samples. 

The superconducting ground state breaks both color and electromagnetic gauge symme- 
tries. It would therefore seem that excitations above the condensate can only be character- 
ized by the unbroken Q charge. At first sight, this would imply that electrons and unpaired 
quarks carry only Q charge and must therefore be assigned only a fig chemical potential. If 
this were indeed the case, it would be impossible to neutralize the 2SC phase in the bulk. 
This is because the condensate is Q neutral, but has color and electric charge that cannot 
be neutralized by particles with only Q charge. The resolution to this puzzle lies in noting 
that our expectation to assign only those charges that are unbroken by condensation to 
excitations applies only to excitations above a charge-neutral ground state. In this case, 
charges associated with broken gauge symmetries are easily delocalized and transported to 
the surface by the condensate. It is however important to note that only the excess bro- 
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ken charge resides on the surface. In describing particles that make up the charge neutral 
ground state we must use vacuum quantum numbers. In this case, the individual charges 
are localized on the particles in the bulk. Therefore, in what follows we treat electrons and 
unpaired quarks as carrying their vacuum charges in our description of the neutral ground 
state. 

Thermodynamics: We begin with the NJL Lagrangian jTU], |TT], |T2|, [TB|, |T^ supplemented 
by both a diquark interaction and the t'Hooft six-fermion interaction which reproduces the 



anomalous Ua(1) symmetry breaking present in QCD [15|. Explicitly, 
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a=0 '- 

-Gd [detij (1 + 275) Qjp + deti^ (1 - ^75) Qjp] Sap 



C 



^15^1' j'k£a'/3'"/<lj'/3' 



(2) 



where rriij is the diagonal current quark matrix, and the spinor q'-^ = Cq'^, where C is the 
Dirac charge conjugation matrix. We use 7 for color (r= red, b= blue, and g= green) 
indices, and i,j,k for flavor {u= up, d= down, and s= strange) indices throughout. The 
chemical potential matrix is diagonal in flavor and color, and is given by 



(3) 



where /i is the quark number chemical potential. Since the couplings Gs, Gd, and Gdiq 
are dimensionful, we impose an ultra-violet three-momentum cutoff. A, and results are 
considered meaningful only if the quark Fermi momenta are well below this cutoff. The 
values of the couplings are fixed by reproducing the experimental vacuum values of /^r, m^r. 



rriK, and m^/ as in Ref. |T^. For the most part, we discuss results obtained using 



GsA' 



niQd = 5 MeV , 



140 MeV 



A = 600 MeV 



1.84, GdA^ = 12.4^ 



and Gdiq = ^Gs/^. 



(4) 



In vacuum, the effective four-fermion interactions in the qq and qq channels are related by 
a Fierz transformation, hence the choice of Gdiq = 3G5/4. Note that Eq. (0) does not 
include the possible presence of a six-fermion interaction due to diquark {{qq)) condensates; 
such interactions have been assumed to result only in a renormalization of the four-fermion 



diquark interaction. In the mean field approximation, the thermodynamic potential per unit 
volume is given by 

I , 2 

I A 



k 1 ^Gdiq 



i=u,d,s 
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(5) 



where {qiqi) {i = u, d, s) is the quark condensate, and the term Qq ensures that the zero 
density pressure, P = —Q, of non-superconducting matter is zero: 



Qo = 2Gs E (9.9.)o - "^Go {uu), (dd)^ (ss), + 2N,Y.J T^lyA'+P^ 



,2 I „2 



(6) 



i=u,d,s I ^ 

where < g^gj >o denotes the value of the quark condensate at zero density. The gap matrix 



A'^^ = 2Gdiq Uatl^e'^'e^^^^qf^ 



(7) 



features three non- vanishing elements. Using the standard notation of denoting A'^'^ through 
the flavor indices i and j, we have 



A,, = A-, 



A„. = A-^^ and A„, = A'" . 



(8) 



This corresponds to the ansatz in Ref . , except that color sextet gaps (symmetric in both 
color and flavor) are ignored. Inclusion of the sextet gaps modifies our results only slightly, 
because such gaps are small 0. Note, however, that we have removed the degeneracy 
between A^^ and A^^^ in order to explore phases in which these gaps may not be equal. 

The quasiparticle energies A, may be obtained by diagonalizing the inverse propagator. 
Equivalently, Aj are the eigenvalues of the (72 x 72) matrix 



D 



-7O7 ■ p - Mi70 + AZ7075C 

70«75A -7^7^ ■ p + Ma"" - fiic 

where Mj are the dynamically generated quark masses and A is given by 



(9) 



Equations (H) through (^) enable a consistent model calculation of the thermodynamics of 
superconducting quark matter as a function of the chemical potentials /x, /ig, /is, and /ig 
at arbitrary temperatures. For a given set of these chemical potentials, the dynamical (or 



constituent-like) masses Mi and the gaps Ajj are determined by the solutions of equations 
that result from extremizing Vl with respect to {qiqi) and Ajj, respectively. 

The phases are labeled according to which of the three gaps are non-zero: (1) Normal 
phase: all gaps zero, (2) 2SC phase: only A^^ is non-zero, and (3) CFL phase: all gaps are 
non-zero. Where needed, we add electrons simply by noting that /ig = — /ig, and include 
their free Fermi gas contribution to Vl. 

Recall that the ground state is Q neutral, i.e., riQ = —dQ/dfiQ = 0, which is a consequence 
of the fact that the condensates are Q neutral. Quasiparticles carrying Q charge are massive 
with m ~ A . In addition, the Q susceptibility Xq = dnq/djiQ ~ 0; in fact, the free energy 
density is independent of /ig at zero temperature. This is because to generate Q charge in 
the ground state, we must break a pair and the energy cost is of 0[A]. In contrast, the free 
energy density depends on fiQ,fi3, and fis, and, the corresponding individual susceptibilities 
do not vanish. For a physical rris of order 100 MeV, there is no apriori reason to expect 
equal numbers of u, d, and s quarks in the CFL phase. The pairing ansatz in Eq. (Wj and 



the arguments of Rajagopal and Wilczek |T^, however, guarantee that 



nrd = , Ubd = Ugs , and rirs = rib^ , (10) 

or equivalently, that 

nu = nr, nd = ng, and Us = rih , (11) 

where Uai is the number density of quarks with color a and flavor i, and Ua {ui) is the net 
number density of color a (flavor i) . Pairing by itself does not enforce either color or electric 
charge neutrality. The strange quark mass induces both color and electric charge in the CFL 
phase. We are, however, at liberty to adjust the chemical potentials /is and /is to enforce 
color neutrality. Moreover, since the pairing ansatz enforces Q neutrality, enforcing color 
neutrality automatically enforces electric charge neutrality at fig = 0. In contrast, the 2SC 
phase requires a finite /ig to satisfy electric charge neutrality and hence admits electrons. 

We turn now to discuss results, beginning with those at temperature T = 0. In Figure |l|, 
we show the dynamically generated or constituent d and s quark masses Mj (left panel) and 
the pairing gaps Aij (right panel) as functions of fi in the CFL and 2SC phases. The u 
quark mass, which tracks the trend of the d quark, is not shown for the sake of clarity. The 
dark (light) curves refer to the case in which color and electric charge neutrality is (is not) 




FIG. 1: Dynamically generated masses and pairing gaps in the CFL and 2SC phases at zero 
temperature from NJL model calculations. Dark (light) curves refer to results when color and 
electric charge (C &: Q) neutrality is (is not) imposed. 

imposed. All masses decrease with increasing /i, since all of the {qq) condensates decrease 
with 11. Note that the requirement of color and charge neutrality has a larger effect on the s 
quark mass in the 2SC phase than in the CFL phase. This is because neutrality in the 2SC 
phase requires a large and negative electric charge chemical potential. In the discussion that 
follows, we will show that /ig ~ — m^/2/i in the 2SC phase. Further, since fig = fi — fig/S, a 
large and negative /ig enhances the strange quark density which in turn suppresses the (ss) 
condensate. 

The right panel of Figure shows the various gaps in the CFL and 2SC phases. Imposing 
color neutrality reduces A^^, since the numbers of red and green quarks (equivalently of u 
and d quarks) are reduced relative to the colored case (see the analytical analysis below). 
For the same reason, color neutrality increases the gaps involving strange quarks. These 
trends are broken only when n begins to approach the ultra-violet cutoff A. The strong 
increase of the A^^ gap as /x decreases for matter in which neutrality is not imposed is due 
to the strong decrease in the gaps involving strange quarks. 

In Figure ||, we show the chemical potentials /zg and fj,Q (as functions of fj,) required to 
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FIG. 2: Chemical potentials /ig and fiQ that ensure color and electric charge neutrality in the CFL 
and 2SC phases as functions of the quark number chemical potential fx at temperature T = 0. 
Solid (dashed) curves refer to results of the NJL (simplified) model. 

achieve color and electric charge neutrality in the CFL and 2SC phases. The solid curves 
refer to results of the NJL model calculations. The left panel of Figure |^ shows the pressure 
P versus /i at T = 0. Here the dark (light) curves refer to the case in which color and electric 
charge neutrality is (is not) imposed. Note that the pressure of the color and electrically 
neutral normal phase falls below that of the 2SC phase for all /is shown. The pressure 
differences AP or the free energy density cost necessary to ensure color and electric charge 
neutrality in the CFL and 2SC phases are shown in the right panel of Figure |^. Here also 
the solid curves refer to results of the NJL model calculations. 

In order to gain a qualitative understanding of the results in Figures ^ and ^, we undertake 
an analytical analysis of a simpler model also considered in Ref. 0. In this analysis, we 
consider u and d quarks as massless, and include corrections due to the s quark mass 
at leading order as a shift in its chemical potential. This does not properly account for 
the shift in energy due the strange quark mass for states far away from the Fermi surface, 
but is a consistent approximation in this context since we are primarily interested in the 
leading order cost of enforcing neutrality. Further, we assume that all gaps, including those 



8 




FIG. 3: Left panel: The pressure P versus quark number chemical potential /x in the CFL, 2SC, 
and normal phases at temperature T = 0. Dark (light) curves refer to the case in which color and 
electric charge (C &: Q) neutrality is (is not) imposed. Right panel: Pressure differences AP or 
the free energy density cost required to ensure C & Q neutrality in the CFL and 2SC phases at 
r = 0. Solid (dotted) curves are results of the NJL (simplified) model calculations. 

involving the s quark, are independent of ^, and that both the gaps and are weak 
functions of the chemical potentials. With these assumptions, and to leading order in A 

^CFL — ^rgb ~t~ ^^rg ~t~ ^rs ~t~ ^gs 5 

^^rb = (/irb + 3AVrb) ' ^nd figb = (^gb + 3AVgb) , (12) 

where we have written the free energy of the CFL phase in terms of the 3x3 block involving 
ru — gd — hs quarks, and three 2x2 blocks involving rd — gu, rs — hu and gs — hd quarks, 
respectively. Each of the three 2x2 blocks is rigid in the sense that the free energy is 
unaffected by differences in chemical potentials of quarks that pair in a given block [p!6| . 
The free energy depends only on the average chemical potential. The 3x3 block does 

^ Corrections to A due to m^. arise at 0[m1/ \ for a detailed discussion, see Ref. Jl7[ . 
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not exhibit this rigidity; here, the quasi-particle energies depend on the sphtting between 
the chemical potential characterizing the ru, gd and bs quarks. Chemical potentials that 
characterize the free energy of the 3x3 block are given by 

/.i=/i + ^, /.. = /.--^--, and /^3 = /.--^, (13) 

and the common chemical potentials that appear in the free energy expressions of the 2x2 
blocks are given by 

m? VP? 

We note that, in general, pairing between particles with dissimilar masses does not require 
a common chemical potential. Maximal BCS-like pairing requires that the distribution of 
the pairing partners be identical in momentum space. Since we treat the u and d quarks 
as massless particles and account for the effects of the s quark mass through a shift in the 
chemical potential in our analytic analysis, a common chemical potential within each pairing 
block ensures that the aforementioned pairing criterion is satisfied. 

In the CFL phase, the stress induced by the strange quark mass generates color charges. 
In the limit of nearly equal and vanishing light quark masses, the CFL scheme in Eq. (|1T]) 
indicates that we will require only a non-zero fig to achieve color neutrality. This justifies 
why we neglect fiq and fi^ in Eqs. (|T^) and (|T^. To leading order in the parameter ml/fi, 



and assuming that the differences between the various gaps are small and /i, ms-independent, 
we find that 



2 

mz 



mt mlA'^ 



(15) 



by requiring dflcFL/d^s = 0. Since n„ = and hence Ur = Ug when there are no electrons, 
fis^CFL) = identically in the CFL phase at zero temperature. Naively, Eq. ([T5|) would 
imply that the free energy density cost of enforcing color neutrality in the CFL phase is of 
0[ml^^]. However, we find that such contributions are absent due to cancellations. This 
result (see the lower-most dotted curve in Figure Q), with ms and A of the NJL model as 
inputs, provides an excellent approximation to the exact NJL result. Utilizing Eq. (0), we 
find an analytic estimate for the free energy density cost in the CFL phase: 

m^ m^A^ 



9 ' 9 



(16) 
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The lower dotted curve in the right panel of Figure ^ shows that this result is in quantitative 
agreement with the NJL model calculation. 

In the 2SC phase, the pairing phenomenon itself gives rise to color charges. The 2SC 
thermodynamic potential, to leading order in the gap and consistent with the approximation 



scheme described earlier, is 



rugd 
^free 



^rugd ~l~ ^free ; 

1 / 4 I OA 2 2 

TT^ (/^bu + A'-bd + A'-bs + l^vs + f^gs + f^Q 



(17) 



The chemical potentials appearing above are defined by 



/^rugd 



/^bs 



A^bu 



2/i8 , 



3 



2fi ' 



and 



/i]-S /^gs 



_ 2/i8 HQ 

A^Q m'^s 
VS 3 2/i ■ 



;i8) 



The condition to ensure color neutrality, dQ2SC / dfis 

1 A2 



+ o 



0, yields 



1^' 



(19) 



3^3 /i 

where we have used a common value of A (independent of /i) in the analytical analysis. 
Note that fis{2SC) does not depend on at leading order. Since pairing in the 2SC phase 
involves red and green quarks, it does not induce a color 3-charge; hence /is (25(7) = 0. 
However, electric charge neutrality in the 2SC phase requires an adjustment due to the 
magnitude of m^. At leading order in a expansion, we find 



IIq{2SC) 



1 mi 



+ O 



ml A^ml 



(20) 



2/i 3 /i 

by setting dQ2Sc/dfiQ = 0. As in the CFL phase, the free energy density cost of enforcing 
color neutrality in the 2SC phase is small, because 0[A'^fi'^] terms cancel and the free energy 
density begins to change at 0[A^]. Similarly, we find that there is no cost for enforcing 
electric charge neutrality in the 2SC phase at 0[/t^]. Using the results in Eqs. ([T9|) and (pO|), 
the free energy density cost of enforcing color and electric charge neutrality becomes 

^^2SC = ^2Sc(/i8, fJ'Q) - ^2SC{0, 0) 



' 4 4A2m2 4 An ^ 



m! A^ A^m\ 



9 ' 9 
/t^ /t"^ 



/i^ 



(21) 
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Although the free energy density costs of enforcing color neutrality in the CFL and 2SC 
phases are of the same order, the cost in the 2SC phase is numerically larger. This is in part 
due to the larger strange quark mass in the 2SC phase and because the free energy density 
cost due to A and dependent terms add in the 2SC phase. The analytical results in 
Eqs. (|T^), (pn|), and (PT|), shown as dotted curves in Figures]^ and 0, compare well with the 
results of the NJL model. 

Phase diagram at finite temperature and lepton content: In the proto- neutron star 
context, matter is subject to stresses induced by finite temperature and lepton number 
chemical potentials |jl8|. Since electrons have both electric and lepton number charges, 



fie = -/iQ + fJ^Le , (22) 

where /iLe = f^ue is the chemical potential for electron lepton number. In order to explore the 
effects of a finite neutrino chemical potential and finite temperature on the superconducting 
phases, we employ the NJL model, Eq. (13) with extensions to include neutrinos and electrons. 
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FIG. 4: Cross-sectional views of the T — fi — fiy^ phase diagram at the indicated values of /x and 
/i^e- In both panels, the dark curves show the phase boundaries, while the dotted curves show 
contours of constant baryon density. 

Figure || shows representative cross-sectional views of the T — /x — fiy^ phase diagram. 
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The left panel displays results at fixed /x = 460 MeV. With increasing temperature, a 
first order transition occurs from the CFL phase to the 2SC phase with ud pairing. For 
A^Le = 0, the transition occurs at T ~ 17 MeV. The corresponding baryon densities are 
n-B = {nu + nd + ns)/3 = 1.06 fm~^ in the CFL phase and ub = 0.94 fm~^ in the 2SC phase. 
At zero neutrino chemical potential, an analytic estimate of the critical temperature for 
the CFL-2SC transition can be obtained by assuming that the gaps in the 2SC and CFL 
phases are nearly equal to each other and to their zero temperature values. We find that 



TT 



4 (5A 2 9 6ms 
^Bcs + ^^^ 



,BCS 



with 



where Aq is the zero temperature gap, 6 A and 6ms are the differences between the gaps and 
the strange quark masses in the 2SC and CFL phases, respectively. Tbcs is the temperature 



at which the gap in the CFL phase would vanish, assuming that A(T) = Apy 1 — (T/Tbcs)^- 
Note that, in general, gaps involving strange quarks do not vanish at the transition, i.e, the 
phase transition is first order. This is because at leading order, the critical temperature 
Tc ~ Ao/2.22 is less than that for the second order BCS transition, Tbcs — Ao/1.76. It 
is clear from Eq. ( ^3] ) that contributions to due to 6 A and 6ms can easily alter this, 
allowing for a BCS like second order transition. If the magnitude of the gap in the 2SC 
phase is larger than that in the CFL phase, is lowered and the transition becomes more 
strongly first order. 

Accommodating a finite lepton number in the CFL phase is expensive, because the re- 
quirement of color and electric charge neutrality in this phase excludes electrons. At T = 0, 
the transition from the CFL to 2SC phase occurs ai file — 150 MeV. The latent heat density, 
TA{dP/dT), lies in the range (2-15) MeV/fm^ along the boundary of the first order phase 
transition. With increasing temperature, the critical lepton chemical potential at which the 
CFL-2SC transition occurs decreases. This is because the gaps in the CFL phase decrease 
with increasing T; consequently, unlocking occurs at smaller fi^,^. 

The right panel in Figure ^ shows the phase boundaries at fixed fi^^ = 200 MeV. For this 
neutrino chemical potential, the CFL phase is preferred above /i = 460 MeV. For low (high) 
values of fi, the region of the CFL phase shrinks (expands) progressively to lower (higher) 
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values of T and jd^,^. In contrast, the 2SC-Normal phase boundary is relatively unaffected 
by increasing values of yU (in the range relevant for proto- neutron star studies), although 
minor variations do occur. Note that with increasing temperature, the phase transition 
switches from first to second order. This switch is due to the fact that the A^^ and A^^ gaps 
decrease along the first order phase transition line. When these gaps vanish (at T ~ 25 MeV 
for /ij^^ = 200 MeV), the phase transition becomes second order. In Figure ^, contours of 
constant baryon density are shown by the dotted curves in both panels. Notice that, for the 
values of /i and chosen for display, the 2SC phase supports lower baryon densities than the 
CFL phase. Across the first order phase transition, the density contours are discontinuous. 
We wish to add that the phase in which A^^ 7^ A^^ was found to be thermodynamically 
disfavored in the range of T, /x, and explored here. 

The consequences of requiring local color neutrality in superconducting quark matter 
with and without neutrinos at both zero and finite temperatures are the principal findings 
of this work. Quantitative results, especially those for quark number chemical potentials 
approaching the ultra-violet cut-off in the NJL model used, should be viewed with some 
caution. Notwithstanding this, the basic qualitative features concerning the phase transi- 
tions appear to be generic, insofar as similar trends are found in our analytic analysis that 
employed a simplified model without a cutoff. We also wish to emphasize that the phase 
diagram in Figure ^ requires revision at low values of /i (or low baryon densities) for which 
a hadronic phase is more likely to be favored. 

Kaon condensation: In neutrino-free matter, Bedaque and Schafer have shown that the 
strange quark mass induces a stress which can result in the condensation of neutral kaons 
in the superconducting quark phase |jT9|, Kaon condensation occurs when the stress 
induced by the strange quark mass m1/{2fi) > m^^o, where m/^o is the mass of the neutral 
kaon. In general, the masses of all pseudo-Goldstone bosons receive contributions both from 
the diquark and quark-antiquark condensates. For example, the mass of the neutral kaon is 
given by 0, H, ma 

mj^o = a rriuimd + rris) + xi^d + , (24) 

where a ~ A^/yU^ and x ~ (qq) ■ At asymptotically high density, where the axial U{1) sym- 
metry is restored and the quark-antiquark condensate vanishes, the dominant contribution 
to the masses is from the diquark condensate. In this case, the kaon masses become small, 
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of order 10 MeV, and kaon condensation is robust. At the densities of relevance to neutron 
stars, the situation is less clear because a finite {qq) can potentially result in larger masses 
for the kaons and disfavor meson condensation |]22|, . Pending a detailed investigation of 
this question within the NJL model, we assess here how kaon condensation can affect the 
structure of the color-neutral phase by assuming that the kaon mass is small compared to 
m^/(2/i). This corresponds to near maximal kaon condensation. In this case, the meson 
contribution to the pressure is easily computed and is given by ||2^ 



Pro = \fl 



ml 



1 + 



(25) 



2/i 

where /^r ~ /i is the pion decay constant in the effective theory describing Goldstone bosons 
21|. The leading contribution to the pressure from the Goldstone bosons is of order m^. 



Including this contribution we find that at zero temperature and neutrino chemical potential, 
the phase transition between the 2SC and CFLKO phases occurs at a value of /x which is 
lower by 16 MeV compared to the case without kaons. 

At finite neutrino chemical potential, CFL quark matter contains a condensate which 
admits electrons, even at zero temperature In the presence of electrons, charged kaon 
condensation lowers the free energy density cost for accommodating lepton number. In- 
cluding the electron and contributions to the pressure {Pk+ = fllAc+l'^)^ ^^"^ solving 
self-consistently for the condition of charge neutrality, the pressure is increased by about 
2 MeV at = 200 MeV. These results indicate that the extent of the meson condensed 
CFL phase in the phase diagram is likely to be enlarged, but only by a few percent. In 
our analysis, the order of the phase transitions between the CFL and the 2SC phase was 
not affected by meson condensation. We wish to note, however, that our analysis neglects 
the effect of the meson condensate on the quasi-particles themselves. This feedback, which 
can alter the quark contribution to the free energy, must be explored before quantitative 
conclusions regarding the role of meson condensation on the phase diagram can be drawn. 
This warrants further investigation and is beyond the scope of this article. 
Astrophysical Implications: The T— /i— /ij,^ phase diagram offers clues about the possible 
phases encountered by a neutron star from its birth as a proto-neutron star (in which 
neutrinos are trapped) in the wake of a supernova explosion to its neutrino-poor catalyzed 
state with ages ranging from hundreds of thousands to million years. In earlier work, some 
aspects of how a phase transition from the normal to the 2SC phase would influence neutrino 
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transport in a newly-born neutron star were explored To date, detailed calculations 
of the evolution of a proto-neutron star with quarks have been performed for the case in 
which only the normal phase was considered Our findings in this work indicate that 
the core of a proto-neutron star may well encounter a 2SC phase first when matter is hot 
and neutrino-rich before passing over to a CFL phase. 

Conclusions: Color and electric charge neutrality in the superconducting quark phases re- 
quires the introduction of chemical potentials for color and electric charge. The magnitudes 
of these chemical potentials are sub-leading in fi. The corresponding free energy density 
costs are small and independent of fi at leading order with the free energy density cost for 
neutrality in the 2SC phase being significantly larger than that in the CFL phase. Conse- 
quently, and in agreement with Ref. 0, we find that the bulk 2SC phase is less likely to 
occur in compact stars at T = and /i^^ = 0. In the NJL model, a small 2SC window 
does exist at relatively low baryon density. However, since this window occurs at very low 
density it is likely to be shut by the hadronic phase. If homogeneous quark matter were 
to occur in neutron stars, it seems likely that with increasing fi a sharp interface would 



separate hadronic matter and CFL quark matter [27]. We note, however, that we have only 
considered homogeneous phases in this study and it is possible that less symmetric hetero- 
geneous phases may well be favored for chemical potentials of relevance to neutron stars. 
Examples include the CFL-Hadron mixed phase and crystalline superconductivity p5|. 
These possibilities alleviate the cost of enforcing color neutrality, since in these cases it is 
only a global constraint. In such phases, however, energy costs associated with gradients in 
particle densities must be met. 

The value of the diquark coupling employed in this work predicts gaps on the order of 
100 MeV at /i ~ 500 MeV. The relationship between the (qq) and (qq) condensates in 
medium obtained by employing the mean-field gap equations with the couplings set by the 
Fierz transformation in vacuum may differ from that obtained in a more exact treatment 
of the NJL model. Lacking experimental guidance on their values in medium, we have 
studied the influence of moderate changes to the diquark coupling (but within the mean- 
field approximation) on the predicted phase structure. For example, using Gdiq — Crg, we 
find that gaps are increased by about 20% in both the CFL and 2SC phases relative to the 
case with Gdiq = ^Gs/'^ predicted by the Fierz transformation in vacuum. The extent of 
the low density region in which a 2SC phase is favored over the CFL phase is not greatly 
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changed and the 2SC phase continues to be favored over the normal phase. 

Had we ignored the differences in the strange quark mass between the charge-neutral 
normal phase and the charge-neutral 2SC phase, the normal phase would be favored over 
the 2SC phase at low density when A < m^/4yU However, we find that the chemical 
potential for the strange quarks is larger in the charge-neutral 2SC phase than that in the 
normal phase. Consequently, (ss) is reduced and the lighter strange quarks in the 2SC phase 
contribute more to the pressure. It is this feedback that tips the balance in favor of the 2SC 
phase when A < m^/4/i. However, we do not expect this trend to remain intact for larger 
variations in the couplings. Obviously, it is always possible to reduce Gdjq and increase Gs 
so as to allow for the existence of a normal phase at lower density. At finite temperature 
and neutrino chemical potential, the CFL phase becomes less favored both because of its 
small specific heat and because of its exponentially suppressed (by the factor exp(— A/T)) 
electron number density, which makes the free energy density cost of accommodating lepton 
number large. In contrast, the 2SC phase has a larger specific heat and easily accommodates 
electron number, and is therefore the favored phase at finite temperature and lepton number. 

The inclusion of Goldstone bosons in the CFL phase tends to extend the region in the 
T—fJ'Ue plane where the CFL phase is favored, since Goldstone bosons contribute significantly 
to the specific heat and also allow for the presence of electrons. In the absence of Goldstone 
bosons, a first order unlocking transition occurs from the CFL phase to the less symmetric 
2SC phase with increasing lepton chemical potential. When the temperature is sufficiently 
high, the phase transition switches from first to second order due to the fact that the A^^ 
and Ads gaps decrease along the first order phase transition line and eventually vanish. 

As discussed above, different phases of color superconducting quark matter are likely to 
be traversed by the inner core of a proto-neutron star during its early thermal evolution. The 
task ahead is to study how these phases and transitions between them influence observable 
aspects of core collapse supernova, neutron star structure, and thermal evolution. 
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